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Complete classification of purely magnetic, non-rotating 
and non-accelerating perfect fluids. 
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Faculty of Applied Sciences TW16, Gent University, Galglaan 2, 9000 Gent, Belgium 

(Dated: February 7, 2008) 

Recently the class of purely magnetic non-rotating dust spacetimes has been shown to be empty 
(Wylleman, Class. Quant. Grav. 23, 2727). It turns out that purely magnetic rotating dust models 
are subject to severe integrability conditions as well. One of the consequences of the present paper 
is that also rotating dust cannot be purely magnetic when it is of Petrov type D or when it has 
a vanishing spatial gradient of the energy density. For purely magnetic and non-rotating perfect 
fluids on the other hand, which have been fully classified earlier for Petrov type D (Lozanovski, 
Class. Quant. Grav. 19, 6377), the fluid is shown to be non-accelerating if and only if the spatial 
density gradient vanishes. Under these conditions, a new and algebraically general solution is found, 
which is unique up to a constant rescaling, which is spatially homogeneous of Bianchi type VIo, has 
degenerate shear and is of Petrov type I(A/°°) in the extended Arianrhod-Mclntosh classification. 

The metric and the equation of state are explicitly constructed and properties of the model are 
briefly discussed. We finally situate it within the class of normal geodesic flows with degenerate 
shear tensor. 



PACS numbers: 04.20.-q, 04.20.Jb, 04.40.Nr 

1. INTRODUCTION 

For a given spacetime geometry, the electric and mag- 
netic parts of the Weyl tensor Caked w.r.t. some unit time- 
like congruence u" are pointwise defined by 



Eab = Cacbdu'^U , 
1 



Hn 



- 7nny^ c a 



(1) 

(2) 



Vabcd being the spacetime permutation tensor. Eab 
and Hab are traceless and symmetric tensors satisfying 
Habu'' = EahU^ = 0, and determine the Weyl tensor 
completely [H y, y|. They were first introduced (for 
the vacuum Riemann tensor) by Matte [jj] when search- 
ing for gravitational quantities playing an analogous role 
to the electric and magnetic field in classical electro- 
magnetism. Using the decomposition (Q, the Bianchi 
identities take a form analogous to Maxwell's equations 
for the electromagnetic field ''^. A non-conformally flat 
spacetime for which Eab, resp. Hab, vanish identically 
w.r.t. some mq" has therefore been called purely mag- 
netic (PM), resp. purely electric (PE), and its Weyl ten- 
sor is said to be PM, resp. PE, w.r.t. uq°-. As Eab and 
Hab (w.r.t. any u°') are diagonalizable tensors, the Petrov 
type of PE or PM spacetimes is necessarily I or D, and 
in each point mq" is a Weyl principal vector (which is 
essentially unique for Petrov type I, and which is an ar- 
bitrary timelike vector in the plane of repeated principal 
null directions for Petrov type D) |^ |y|. The PM and 
PE property may be characterized independently of u°, 
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as follows. Defining 



= Ea 



iHab, (3) 

the quadratic, cubic and 0-dimensional invariants /, J 
and M of the Weyl tensor [2, i8|| can be written as [Ij, i^] 



I = Q\Q\ 



- li + li + ll = -Hhh + hh + hh) 



■2 I 72 

{E\E\ - H\H\) + 2iE\H\, 



J - Q\Q\Q\ 



= il + i"^ 



2 ^ '3 

lb 



I3 — 0I1I2I3 



M=--6 



Tpa rpb rpc o rpa rrb tjc 

—i{H'^i,H cH'^a ~ 3E"'i,E cH'^a), 






(4) 



(5) 



(6) 
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where the solutions li ofl^ — ^lI—^J — are the non-zero 
eigenvalues of Qab for any u". Hence, a spacetime is PE 
(PM) if and only if all k are real (imaginary), or ^ if and 
only if M is real non-negative or infinite and I is real pos- 
itive (negative). For PE or PM spacetimes, M = if and 
only if the Petrov type is D; the respective types in the 
algebraically general case were logically denoted I{M+) 
and I{M°°) in the extended Petrov classification of Ari- 
anrhod and Mcintosh [1^, where I{M°°) corresponds to 
J = 0, i.e. to one of the Z, being identically zero. 

Whereas large and physically important classes of 
examples exist for PE spacetimes (for example all the 
static spacetimes are purely electric |47j|), only a few 
examples exist of the purely magnetic ones. This is 
particularly true for the vacuum solutions (with or 



without A term), where no PM solutions are known at 
all. This has lead to the conjecture that PM vacua do 
not exist Q, but so far this has only been proved for 
Petrov types D and /(M°°)0| or when the timelike 
congruence u" is shear- free [T^j non-rotating [13, [lJ|, 
geodesic [T^ , or satisfies certain technical generalizations 
of these conditions [la, [l3 ■ In [iM uM ^be non-existence 
of shear-free or non-rotating PM models was generalized 
to spacetimes with a vanishing Cotton tensor. As a pos- 
itive example on the other hand, the metric constructed 
in [231 turns out to be a PM kinematic counterpart 
to the Godel metric [2l|, but its source is unphysical 
as the Ricci tensor is of Segre type [11, ZZ]. In 22], 
PM locally rotationally symmetric (LRS) spacetimes 
were shown to belong to either class I or III of the 
Stewart-Ellis classification 23], the possible Segre- types 
were determined and the most general metric forms 
were found, exhibiting one arbitrary function and three 
parameters. 

In a cosmological context, perfect fluid models are 
studied, the metric gat being a solution of the Einstein 
field equation with perfect fluid source term 



Rab — -^Rgab 



(/^ -I- p)Ua Ub + PSab- 



(7) 



Here Rab is the Ricci tensor, which is assumed not to 
be of A-type {fi + p y^ 0). In this case, the average 4- 
velocity field m" = Upf° of the fluid plays the role of a 
geometrically preferred unit timelike vector field on the 
spacetime. A possible cosmological constant A has been 
absorbed in the fluid's energy density and pressure by 
means of re-defining /i = /i' -f A and p = p' — A. The 
electric and magnetic parts Q) w.r.t. Wpf" then represent 
the locally free gravitational field (not pointwise deter- 
mined by the matter content via the Einstein equations 
0), while the vanishing of their spatial divergence to- 
gether with the non- vanishing of the so called 'curl' and 
'distortion' parts of their spatial derivatives form neces- 
sary conditions for the existence of gravitational wave 
perturbations of the homogeneous and isotropic FRW 
spacetime |2J, |2^ |2^ ^31 ■ Whereas the electric part 
Eab is the general relativistic generalization of the tidal 
tensor in Newtonian theory [23], the magnetic part Hab 
has no Newtonian analogue and it's role is still poorly 
understood (see however |29l ISOL l81| . 

This motivates the attempt to construct purely mag- 
netic perfect fluid models with uq"" ~ Upf° in the above. 
A solution to {Tj) of this kind will then be called a PMpf. 
The only known perfect fluid which is purely magnetic in 
a certain region of spacetime, but which is not a PMpf 
in the above sense (uo° 7^ Upf°), is the Lozanovoski- 
Mclntosh metric [33 ]. 

The present study is restricted to PMpf models, with 
the exception of section 4.4. First however we conclude 
this introduction with an overview of all performed stud- 
ies and known examples of PMpf's, to the best of our 
knowledge. 



In the case where the Petrov type is D, an exhaustive 
classification of PMpf's has by now been obtained, 
essentially due to two results. Firstly, it was shown most 
recently "SS] that Petrov type D PMpf's are necessarily 
LRS class I or III, whereas secondly the most general 
metric forms [Sg for each of these classes were derived in 
[2^ . Four systematic studies historically precede these 
results: on the one hand LRS space-times in general 
were reinvestigated in [S^ [3g by two different methods, 
hereby also briefly discussing the purely magnetic case; 
on the other hand shear- free [a], resp. non-rotating [33 
PMpf's of Petrov type D were shown to be LRS class 
I, resp. III. Finally the following three distinct PMpf 
solutions are particular cases of the general LRS cases 
above: the p = i/i Collins-Stewart space-time ^3\ and 
Lozanovski-Aarons metric ^], both LRS class III, and 
the LRS class I rigidly rotating axistationary model 
with circular motion of [34| . The last example was 
found within a study (regardless of the Petrov type) of 
axistationary perfect fluids, which concentrated mostly 
on the PMpf subclass. Concretely, it was shown that 
axistationary PMpf's with circular motion necessarily 
have non- vanishing vorticity and spatial-3 gradient of the 
matter density, and that such spacetimes are LRS class I. 

No algebraically general PMpf's have been found yet, 
although they might be of relevance for cosmological 
modelling. The PMpf subclass of irrotational dust space- 
times consists of 'silent' universe models [2^ [43 and 
was investigated in [2^, but the appearance of chains 
of severe integrability conditions (analogously as for the 
widely studied purely electric subclass in the Petrov type 
I case J£j 44J) made the authors conjecture that this sub- 
class might be empty. This was recently proved in [43] 
for the case of dust, regardless of Petrov type and cos- 
mological constant. 

The present study continues this line of investigation. 
Whereas evidence is provided that the class of PM rotat- 
ing dust must be severely restricted as well, we present 
as a main result a first example of an algebraically gen- 
eral PMpf solution, which is both non-rotating and non- 
accelerating. Up to a constant rescaling, this solution 
is shown to be the unique PMpf with these properties. 
Moreover it is spatially homogeneous with degenerate 
shear tensor and, just as the related LRS class III models 
of 37] , it turns out to satisfy the energy conditions in an 
open subset of spacetime. 

The structure of the paper is as follows. In section 2 
we set up the basic variables and equations for PMpf's, 
in a mixed l-f3 covariant/orthonormal tetrad approach. 
After pointing out why solution families of the PM ro- 
tating dust class should be rather poor in number, the 
subclasses of Petrov type D and of vanishing spatial 3- 
gradient of the energy density are shown to be empty in 
section 3. Section 4 provides a characterization of the 
new metric, discusses its mathematical and physical fea- 
tures and situates it within the broader context of [43. 
We end with a conclusion and discussion in section 5. 



2. BASIC EQUATIONS FOR PMPF'S 

We use units such that SttG = c = 1 , Einstein summa- 
tion convention, the signature (— , +, +, +) for spacetime 
metrics gab, and abstract (Latin) index notation 441| for 
tensorial quantities (except for basis vector fields which 
are written in bold face notation) ; round (square) brack- 
ets around indices denote (anti-)symmetrization. The 
perfect fluid field equations may be rewritten as 



Rab = ■7^{pi- + 'ip)uaUb 



1 



{H-p)hab- 



(8) 



Here hab = Qab + UaUb projects orthogonally to the fluid's 
4- velocity field m", 'spatializing' indices of tensorial quan- 
tities by contraction; if this contraction is the identity 
operation for all indices, the tensor is called spatial. For 
covariant differential operations orthogonal to w°, the 
streamlined notation of 4, 21, 4.5J is the most transpar- 
ent. The covariant spatial derivative D and the associ- 
ated curl and divergence (div) operators, acting on vec- 
tors and 2-tcnsors, are defined as: 

DaS^-''e...f = ha'h^ ■ ■ ■ h\hj ■ ■ ■ hf'\/bSP-%...lC,) 

div V = DaV", curl Va = eabcD'' V^. (10) 

div^a^^'^afc, CmlSab^e.diaD'S^''. (11) 

Here Vc is the covariant derivative associated with the 
Levi-Civita connection and Cabc — Vabcdu"^ is the spatial 
permutation tensor. The kinematics of the perfect fluid 
M°-congruence are then described J43| by its acceleration 
iia, vorticity uja — — ^curlwa and expansion tensor 9ab = 
Di^Ub), which are all spatial. Here and in general, a dot 
denotes covariant ('time') differentiation along u". For 
scalar functions / and vectors w° one has jj,|4^ 

[Daiy = Daf-9a''Dbf-eabcU;''D'f 

+Uaf + Uau''Dbf, (12) 

CnTliDf)a = ^fUa, (13) 

curl(/«), = fcUT\Va + eabcD''fv^. (14) 



The expansion tensor is further decomposed as 9ab 



<^ab 



j6h„h, with a„h = £>, 



{aUb) 



\9K,, = Di 



(ab) 



the 



trace- free shear tensor of the fluid and 9 = 9"- a — div u its 
scalar expansion rate. In general, Si^ab) = ha'^hb Si^cd) ^ 
hScdh'^'^hab is the spatially projected, symmetric and 
trace- free part of Sab, while Vu) = ^a Vb denotes the 
spatial projection of Va 27]. 

In a H-3 covariant 'threading' approach, the tenso- 
rial quantities Ua,i-^a,o'ab,9 and Eab, Hab, IJ-,P are taken 
as the fundamental dynamical fields. One focusses on the 
Ricci-identity for Ua and the Bianchi-identities, wherein 
the Ricci-tensor is substituted for the right hand side of 
(|5|). They can be split in constraint equations (involv- 
ing only spatial derivatives) and propagation equations. 
For a general perfect fluid, the Ricci-equations become. 



our notation: 




2 
div (Ta - ■::Da9 - curlwa + 2eabcUj''u'' = 0, 


(15) 


divw — ii^'uJa = 0, 


(16) 


curl (Tafc + £)(aWb) - Hab + '2U{a^b) = 0, 


(17) 


2 

(^{ab) + 1:0 Gab + (Tc{a(^b)'' +^{a^b) 




-D(aUb) - U(aUb) + Eab = 0, 


(18) 


9+\9'^ + CTabCr"'' - ^UJ^UJa " div U 

O 




-UaU" + -(M + 3p)=0, 


(19) 


2 1 

^(a> + g^^Q - (^abi^ + ^CUrllia = 0. 


(20) 



The equations of conservation of momentum and energy 
(contracted Bianchi-identities) are: 



DaP 



-{jl+p)Ua, 
-{li+p)9. 



(21) 

(22) 



For PMpf's Eab — 0, the remaining Bianchi-identities are 

1 



[<7,H]a-iHabUj'' 

div Hn 



-DaP = 0, 



{p,+p)uJa = 0, 
curl Hab - 7;iP'+ P)'^ab 4 



'2u'^£cd{aHb} 



d 



(23) 
(24) 
0, (25) 



H{ab) + OHab - SadaHb)" + Lu'^e^diaHb) — 0. (26) 

Here [5, T]q = eabcS^'^Td^ is the vector dual to the com- 
mutator of spatial tensors Sab and Tab- Note that in 
general, the term Eab in P8|l couples the evolution of the 
kinematical quantities to that of the Weyl tensor, which 
is no longer the case for PMpf's. The constraints H23() 
and H25|l express div£^a = and Eab — 0, respectively. 

In an orthonormal tetrad approach, a specific or- 
thonormal basis of vector fields B = {eg = u, ea} is 
taken. Here and below, Greek indices run from 1 to 3, 
expressions containing these have to be read modulo 3 
(e.g. Xa+ia-i = Xi2 for a = 3) and daf is written for 
the action of e^ on functions /. For spatial tensorial 
quantities, only the components with Greek indices sur- 
vive; one has in particular hap = 5ap, and eap-y becomes 
the alternating symbol on three indices, where we take 
the convention |57l | £123 = ?7i230 = 1- The commutator 
coefficients ^°'i,(. = —l°'cb and Ricci-rotation coefficients 
V^bc of B are defined by 



[eb, Be] = l'\c^a, ^cGb = r^fccBQ. 



(27) 



As for any rigid frame, the lowered coefficients Tabc — 
gad^'^bc = -Tbac and 7abc = QadYbc are one-one related 
by 



1 



labc 



2ra[bc]i Taftc = 7:ilbac + Icab - labc)- (28) 



Combinations hereof, together with fi,p and the compo- 
nents Ha(3 play the role of basic variables. In the present 
paper, we will use|5^ 



ha = ^a+lOa+l ^^a-lOa-l 
— f^a+la+1 ^ CTa^ia-l 



(29) 



Ca+lQ 



9 


= 


r'^o/s or 




9 


= 


TaOQ = O-aa + o^' 


(30) 


1-1 


= 


x(ra+10Q-l + Ta-lOa+l), 


(31) 


LOa 


= 


;^(ra+10Q-l ^ ^a-lOa+l), 


(32) 


Ua 


— 


TaOO, 


(33) 


n^ 


= 


Tq-Iq+IO, 


(34) 


Qa 


~ 


Ta— 1 a— lo -Laa— 1q— li 




ra 


= 


Ta + l a a+1 ^ a a+1 a + 1: 


(35) 


ria 


= 


Ta+la-la Or 




^aa 


= 


7qq+1q-1 = fia+l + '^Q-l- 


(36) 



In the line of the Cartan one- form formalism, the basic 
equations are the commutator relations (|77|l , the compo- 
nents of the Ricci-identities and of the (second) Bianchi- 
identities ^-^ (see e.g. H^ §7). The Ricci-identities 
may be further split into 
(a) the Jacobi-identities (first Bianchi-idcntities) 



l( d 
Q \a b c 



9[a7^c]+7^[fcc7a]/=0, 



(37) 



with 



l(TC|l and 







(EOl), 




1 2 3J ' ^^^ """ \Q a + 1 a-1^ 

(b) the defining equations (|17|l and (|18() for Hap and Ea/s, 
and 

(c) the components of the Einstein equations (|HJl , where 
the (00) component is Raychaudhuri's equation (|19|l and 
the (Oa) component is H15(l . whereas the (a/3) component 
(which is not covered by the Ricci-identity for Ua) reads 



drr 



{af3) 



ddsT'^a-jc + r'^dcr (q^) - r'-Q/r-' 



/3c 



-(^-p)(5„/3. 



(38) 



Note that (|19() - (|23|l are equalities between spatial tenso- 
rial quantities, which are moreover symmetric and trace- 
free in the case where they are 2-tensors. The non-trivial 
components are readily deduced from (j!!|l- Hll|l . I|35|l - I|3()|) 
and 



3. NO-GO RESULTS FOR PM ROTATING DUST 

The identity ((T^ applied to f — p, together with l(^ 
and idJ, yields 



{fj. + p) curl iia + f^abcD^'p u" 



-2pwa- 



(41) 



Thus non-accelerating perfect fluids are non-rotating or 
dust (i.e. p is constant). From this viewpoint, the class 
of non-rotating dust forms the intersection of both pos- 
sibilities, and in |42| it was shown that PM non-rotating 
dust spacetimes (so-called anti- Newtonian universes |21|) 
do not exist. Natural questions which then arise are (a) 
whether PM rotating dust is allowed, and (b) whether 
non-accelerating and non-rotating PMpf's exist. As said, 
we will mainly concentrate on question (b) , which will be 
dealt with in the next section. Regarding question (a) 
we mention here that, just as in the non-rotating case 
[2l] |. PM rotating dust spacetimes are subject to at least 
two chains of severe integrability conditions, built from 
the consecutive propagations along w° of the constraint 
equations (|23|l and H25|) . A detailed analysis will be given 
elsewhere, but it is clear that only a few distinct fami- 
lies of solutions will be allowed, at most. In this section, 
we explicitly show that the answer to (a) is negative for 
the subcases of Petrov type D and of vanishing spatial 3- 
gradient Da^, of the energy density. The result for Petrov 
type D follows implicitly from .3§|] , but the reasoning pre- 
sented here provides a more transparent and direct way 
to prove it. The result for DafJ, = will be used for 
a further characterization of the new metric in the next 
section. 

In the proofs below, we will use that for iia — Q and 
hi — 0, the diagonal components of (fH^ reduce to the 
single equation 



doh2^^{h2~2e)h2+uol. 



(42) 



Note that for Petrov type D, /i +p = is not allowed an 
se by the result of Q. 

Theorem 3.1 The class of PM type D dust space- 
times is empty. 

Proof. Suppose a PM type D dust spacetime 
{M, gabjUa) exists. Taking an orthonormal eigenframe 
B of Hab, we may assume that H22 = ^33. In this case 
the (22), (33) and (23) components of ^, ^ immedi- 
ately give qi = — ri, (T23 = 0, 7133 = ri22 and hi — (such 
that h^ — — /12), while the (12), and (13) components 
of (HHJ), ^^ together with the (2,3) components of ^^ 



DaVp 


= daV0-VsTl^, 


iea 


,d to 






Va 


= doVa+eaM^^V\ 


(39) 


W2 = --(cri3 + ^2), W3 ^ -{ai2 - 

6^2292 + {^J■+ p)(a'i3 - 3UJ2) = 0, 


-ujs) 


(43) 


Jai>Pj 


= daSf3j - Ss{^Tfj^^, 








Sal3 


— doSal3 + 2e^A-(aO'''S'^) . 


(40) 


GH22r3 + ip + p)(o-i2 + 3^3) = 0. 




(44) 



Now from the (22), (33) and (23) components of H17() and 

,0 1 3] + io 1 2]^^^io 3 ij^io 1 2, 
gets 



from 



(cri2 + LU3)q2 + (cri3 - UJ2)r3 = 0, 
(cri3 - UJ2)q2 - {<Jl2 + (^3)r3 = 0, 



(45) 



such that either ctis — U2 and uyi = —'^z, or 92 = ?'3 = 0; 
in the latter case one has (J\z = 3w2, '^n = — 3^3 by (|44(l . 
In both cases the (22), (33) and (23) components of p8|) 
imply (cj3 — a;2)(<^3 +^1^2) = and uj2'^?, = 0, such that 
(T12 = (713 = CJ2 = '-■^3 = 0. Under these conditions, H23|) 
yields 



TOl 



-IShuji, m2 — TW3 = 0; 



(46) 



while from H2U|I and 119|) we obtain in the variables (|30|l . 



daijJi 



-2^22(^1, 

1 



ao022 = Cjf--(/i + 3p) 



6 



92 

-'22- 



(47) 



However, with 



133) and uji ^ Q ^ ^ 

_ 1 , 



p one deduces 
from [^2, ^3]^^ and ^ that 6I22 = i(/i2 + 6*) = 0. Taking 
two 9o-derivatives hereof and using H47II , (|22|l , p = and 
ji+p 7^ 0, yields 6* = /12 = in contradiction with (|42(I .D. 

Theorem 3.2 A non-accelerating PMpf with van- 
ishing spatial 3-gradient of the energy density is 
non-rotating. 

Proof. Let {M,gab, Ua) be a PMpf with iia = DaP = 
0. As non-accelerating perfect fluids are dust or non- 
rotating, we may assume the Petrov type is I, since for 
Petrov type D the result follows from Theorem 3.1. The 
algebraic vector constraint 1^5)) and its covariant deriva- 
tive along u" yield 



^^c{aH''b)iJ' 






- iHabOj" = 0, 



(48) 
(49) 



Projected onto an orthonormal eigenframe B of Hab, the 
components of these equations form a system of 6 lin- 
ear equations in the variables cri2, fis, o'23, ft.2 and /13, 
parametrized by the Haa and LOa , which can only be con- 
sistent if the determinant of the so called extended system 
matrix vanishes. Writing a; „ = i/a+ia+i— -ffa-ia-i, this 
yields 

^1(^2 + ^3)^2^3 + cyclic terms — 0. 

Hence, as the Petrov type is I, at least two of the spatial 
components uja must be zero, say 072 = ^3 = 0. Herewith 
EHl, 63 imply cri3 = cri2 = 0. Then W2 ^ Wy, = Q 
follows from H26() and (|49|l reduces to 



Taking a further covariant time derivative hereof, one 
finds 



{cTcdCr'"^ - SuJcUj'')Habi^'' 

—abccr'^ H d^a — 0. 



2acdH'^ aab^ 



(50) 



Suppose now uji ^ 0, such that 7^ a; oc ei is an 
eigenvectorfield of Hab- Eliminating (T23 from the first 
components of (|5n|) and (|1S|| (the latter being given by 
f23a;i + 'j-'i(a;3 — X2) — 0), and then eliminating /i2, resp. 
/i3 from the result by means of (|5J) (i.e. h2X3 + h3X2 — 0), 
one finds 

ufx2X3{x2 - X3){Awlxl + hlxf) 

= Ujfx2X3{x2 - X3){4:WJxl + h\x\) = 0. 

Hence X3 = X2 (i.e. Hn = 0), such that hi = (723 = 
by ^ and 148() . while the first component of H26() yields 
Wi — — wi/2. Now the (11) and (1) components of l|^ 
and H24|l respectively read 



("-22 - "-33)2:2 - ^(p + p) = 0, 
(gi +ri)x2 +uji{fi + p) = 0, 



(51) 
(52) 



6aabH'^ujc - (JbcH"^^, 



be, 



0. 



which can only be consistent for X2 7^ if 
wi(»^22 - ^^33) + I (91 + ri)h2 = 0. On the 

other hand, the (23) component of H17|l yields 
wi(n22 - n33) - (gi -|- ri)/i2 = 0, such that 7122 = "33, 
whence /12 = from l|51|l . again in contradiction with 

Remark. This generalizes the result of |Ij| . The rea- 
soning there was similar, but shortcuts were available 
due to the non-existence of Petrov type D j^j and type 
I(M°°) [B vacua. 

From lial and theorem 3.2 it follows: 

Corollary. PM dust spacetimes with vanishing spatial 
3-gradicnt of the energy density do not exist. 



4. ALGEBRAICALLY GENERAL PMPF'S 

Looking at H41|l . (|20|l and (|21|) . one sees that non- 
rotating (non-vacuum) perfect fluids in general obey 
^abcD'' n ii'^ — eabcD^pD'^P = 0, i.e., either (a) the ac- 
celeration Ua and the spatial 3-gradient of the pressure 
DaP are 7^ and proportional to the spatial 3-gradient 
of the energy density Da^i, say DaP = fDaP, / 7^ 0, (b) 
iia = DaP = or (c) Da^i = 0. Hence, a non-rotating per- 
fect fluid exhibits a barotropic equation of state p = p(/i) 
if and only if either (1) the fluid has property (a) with 
/ = |£ = f = -^^ ^ 0, (2) the fluid has the proper- 

ties (b) and (c), or (3) it is dust (constant p, i.e. ^^ = 0). 
For non-rotating PMpf's, it was shown in |43| that the 
dust case (3) is impossible. On the other hand, non- 
rotating type D PMpf's were fully classified in jsjl- Such 



spacetimes were shown to belong to class (2) above and 
turned out to be spatially homogeneous and locally rota- 
tionally symmetric (LRS) of class III in the Stewart-Ellis 
classification. The general metric form as well as the cor- 
responding equation of state were explicitly constructed, 
and the solutions satisfied the energy conditions in an 
open region of spacetime. This result motivates the inves- 
tigation in the present section of the existence of Petrov 
type I non- rotating PMpf's obeying (2). It is proved 
that for non-rotating PMpf's (b) and (c) are actually 
equivalent, and that non- rotating and algebraically gen- 
eral PMpf's exist for which the spatial gradients of mat- 
ter density and pressure vanish. The solution is unique 
up to a constant rescaling of the metric. This provides 
an affirmative answer to question (ii) at the beginning of 
the previous section. The spacetime is characterized in 
Theorem 4.2: it is found to be of Petrov type I{M°°), 
has a degenerate shear tensor and is spatially homoge- 
neous of Bianchi type VIq. Both the metric and equa- 
tion of state are constructed, and the physical behavior 
of the model is briefly discussed. Finally, the solution is 
situated within a larger class of perfect fluid models by 
dropping the purely magnetic condition. 



4.1. Characterization. 

Suppose {M,gab,Ua) is a Petrov type I non-rotating 
PMpf, and project the basic equations w.r.t. an orthonor- 
mal eigenframe B of aab, i-S-i <^a+ia-i = 0. Then the 
(aa) components of (|17|l become algebraic and very com- 
pact in the riQ-variables: 



Inserting (|5^ - l(K5|l into 



Haa — na+lha+l — Ha-lha-l- 



(53) 



On the other hand, one can solve the (Oa) field equa- 
tions and the (a + 1 a — 1) components of (|17|l for 
daha+i,daha-i, giving 

daha+l = ^daO — Ha+la-1 — ?'q^q-1 — ^Qaha+i^A) 
daha-1 = daO - Ha+la-l+'iraha-l+qaha+l^i^b) 

while (|18() and (|19() combine to prescriptions for the evo- 
lution of the variables (|30|l : 



doOaa = -0^^ + daila + il-l + Qa+lUa+l - ra^lUa-l 



-^(M + 3p). 



As shown below and in the appendix, we may take Wa = 
for the cases DaP = and DaV = Wa = under study. 
Then the (aa) components of (j^ reduce to 



+ 7:ha{Ha-la-l — Ha+la + l)- (57) 



a-1 
a a — 1 



and 



a + 1 
a a + 1 



one gets 



^O^Q — 7;{f^a — ha+l — 0){qa — Ua) ~ Ha + la~l, 
9ora = -{ha-1- ha-0){ra+Ua) - Ha + la-l{58) 



while _ , _, _, yields 

\0 a + 1 a-l'-' 
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doUaa = -l^OUaa + ^ C^a-l " /^Q+O'^aa, (59) 

or translated for Ha instead of Haa'- 

doUa = --Oria- -{ha^i-ha+l)na 



We will make use of the following 



(60) 



Lemma 4.1 W.r.t. a shear-eigenframe of a non- 
rotating PMpf: 

(a) at most one of the n^ vanishes; 

(b) at most one of the h^ vanishes; 

(c) if iia = 0, then a/i2 + bh^ = for constant a, b 
((a, 6) 7^ (0,0)) implies /ii/i2^3 = 0, i.e. the shear tensor 
is degenerate. 

Proof, (a) Immediately follows from (|60l) and H53|) . 

(b) This is the statement that for non-rotating PMpf's 
the shear tensor cannot vanish 13], as is immediately 
seen from (|17ll in covariant form (or from H53|) and H54|) - 
(|55|l in tetrad form). 

(c) For iia = 0, the difference of the (a -I- 1 a -|- 1) and 
(a — 1 a — 1) components of (|18() reads 



doha — — -(^Q+i — ha-i + 29)ha- 



(61) 



From this and /ii + /i2 + /13 — 0, one deduces 
do{ah2 + bhj,) + ^{29 + hi - ft,2)(a/i2 + 6/13) = a/ii/12, 
from which the result follows. D 

Theorem 4.2 For algebraically general PMpf's, any 
two of the following three conditions imply the third: 



(56) (i) the fluid congruence is non-rotating; 



(ii) the fluid is non-accelerating (i.e. the spatial gradient 
of the pressure vanishes) ; 

(iii) the fluid has a vanishing spatial gradient of the matter 
density. 

Proof. 

(i), (ii) =^ (Hi). See appendix. 

(ii), (Hi) => (i). This was the content of Theorem 3.2. 



{iii),{i) ^ (ii). When 9q^ — one deduces from 
[do,da]^i (or ^) and (EB that also da9 = 0. With 
Hfj+ifj-i ~ 0, (|^ becomes algebraic: 

{Haa — Ha-la-l)na+la + l + {Haa — Ha+l a+l)na-l a-1 
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-{ha+l - ha-l) = 



(62) 



and consists of two independent equations. Elimination 
oi ^1 + p yields 

{h2 - h3){H22 - H33)ni + cyclic terms =0 (63) 

From (|53|l one deduces 

riaHaa = "-a ("-q+i/Iq+1 " '^q-I^q-i), (64) 

{nin2 + ri2n3 + n:>,ni)ha = 

Ua^lHa-la-l — na+lHa+la+l- (65) 

li X = nin2 + ri2'T-3 + n^ni vanishes, then Y = niHn — 
"-2-^22 = "-3-^33 from (|65|l and OqX — ~2Y = from 
(I64II . which is contradictory to lemma 4.1 (a). Hence 
X 7^ and one can use H65|l to eliminate the ha- Doing 
this for H63|l one obtains (ni + n2 + n^)F = 0, with 



F = niHii{H22 — H33) 
This leaves two cases 



cyclic terms 



(I) F^O: 

Calculating OqF + ^9F and substituting for ha one 
gets another polynomial equation G = in n^ and 

Haa- Substituting Ha+la+l = —Haa ~ Ha-la-1 and 

taking the resultant of F and G w.r.t. H^-i a-i yields 

Haa{ni - n2)(?l2 - 'T-3)("-3 " "l) 



X {nin2 + n2n3 + nani)^ 



0. 



Hence e.g. n^ = 712. 
H,i{H22 - H33){ni - 
G = would read H. 



Inserting this in F = yields 
712) = 0. If Hii were zero then 
= 0, hence 712 = n^, = 0, 



22 "-2 
again in contradiction with lemma 4.1 (a); if rii — n2 

were zero then G = would read n2(-ffii —-^22) (^22 ^ 
H33){H33-Hii) = 0. Thus F = leads to Petrov type 
D (and hence to the models found in |33), contrary to 
our type I assumption. 



(II) m + n2 + ns — 0: 

Applying do to this equation one finds 

(ft,2 - /i3)ni + (/13 - hi)n2 + (hi ~ ft.2)«3 = 

or, by (jn3). 



Pi = n\Hix + nlH2 



nlH33 = 0. 



(66) 



(67) 



Calculating doPi + ^9Pi and again eliminating the ha 
by H65I) one obtains a further polynomial equation P2 
in the Ua and Haa- Substituting iJa+ia+i ~ —Haa — 
Ha-ia-i and eliminating Ha^ia-i from Pi and P2 
now yields 

-ffaQ"-l"-2'^3("-l"2 + ?^2»^3 + ?^3?ll)^ = 0. (68) 



Hence e.g. ni = 0, such that 712+71.3 = tt-h = 0, ft,i = 
from (EHl) and iJn = from ((HTI). From hi = and 
(fK^ . if^ one deduces qi = — ri and (72 = 0,r3 = 0. 
Now doq2 — dofs — yields U2{2h2—9) — 713(2/7,2— 0) = 
by means of (|58|l . 

Suppose that 7x2 and 7x3 are non-zero, such that 2/i2 — 
9 = ^39ii == 0. Then di{2h2 - 9) = yields /i2(ri - 
2gi) = by lElJl and di9 = 0, such that qi ^ n ^ 
by lemma 4.1 (b). Now dori = yields /12W1 = 
by H58|l and hence, again by lemma 4.1 (b), Wi = 0. 
Then however, by (|56|) for a = I, one would have fi + 
3p = 0, leading to daP = iia — 0, contradictory to 
the assumption (7x2,1x3) 7^ (0,0). Hence 7x2 = 7x3 = 0. 
Finally, the (23) component of fl^ becomes 7x1(7^2 — 
7x3) ~ 0, such that also 7x1 = by lemma 4.1 (b), and 
the fluid is non-accelerating. 



4.2. Existence and properties. 

From the proof of theorem 4.2 it follows that for 
non-rotating, algebraically general PMpf's with vanish- 
ing spatial gradients of matter density and pressure, one 
of the eigenvalues of Hab is identically zero, such that 
the Petrov type is I(M°°) in the extended Arianrhod- 
Mclntosh classification. The i?afc-eigenframe {u, 60,} be- 
ing also an eigenframe of the shear tensor, the lat- 
ter is degenerate in the plane perpendicular to the 0- 
eigendirection of Hab- When ei spans this eigendirec- 
tion, i.e. for Hu = hi = (such that iJ33 = — -^22 7^ 
0,/i3 = — /i2 7^ 0), the following equations have by now 
been established: 



daP = dafJ. = da9 = 0, 

Ua — ^a — O'a+la-l — 0, 

Wa = 0, 

ni = 77,2 + ns = (i.e. n°'a 

qi+ri = 92 = r^ = 0, 



while H53I) and H62|l reduce to 



H22 = "2/7,2, 
1 



Ho 
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(m+p). 



(69) 

(70) 

(71) 

7X11 = 0), (72) 

(73) 



(74) 
(75) 



Now from lf75|l and ifS^ . (jS^J one sees that ^2^2 = 
^3/12 = 0, while with ^^, ^^ and ^^ one immedi- 
ately deduces from [dQ,di]9 that also 9i/i2 = 0, whence 
gi = n = by lO and ^- From the {a + la- 1) 

field equations and I _. „ „j, [di,d2]n2 and [93,9i]n2 

one deduces 93 = 7-2 = and 9a?X2 — (alternatively, 
dan2 = daH22 — follows from ifT^ and (|73|l . whence 
^2 = 93 = by (|24|) . but the previous reasoning is more 
generally valid, cf. infra). Finally, the (aa) field equa- 
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tions and H75|) are equivalent to 
1 



2 = T:^22ifiii — ^22) 



M = X^22(^ll + ^22), 
P = ^^22(^11-3^), 



(76) 
(77) 
(78) 



where, in view of the construction of the metric below, 
we re-introduced the variables On and ^22 = ^33- By 
virtue of ((23, (EOJ), ^^ and (|^, the relations (|7^ . 
(|76|l and (|77|l are consistent under propagation along the 
matter flow lines, whereas applying 9o to H78|) yields an 
expression for p in On and ^22^ 



P 



-W- 



22lfii — C722ncil 



){0i 



202 



(79) 



With the above specifications, all basic equations are sat- 
isfied and consistent, implying that corresponding solu- 
tions exist. Since all ^a-derivatives of the commutator 
coefficients vanish the spacetimes will be spatially homo- 
geneous, while (/q = Tq, = and l|72() imply that the 
Bianchi type is VIq (see H^)- 



4.3. Metric, uniqueness and equation of state 

With (|7ni)-(|72I), g„ = r^ = and 6*22 = 6*33, it can be 
easily checked that the normalized vector fields 

eo = u, ei, 

62' = ^(62 +63), 63' = — =(e2 -es) (80) 

are hypersurface-orthogonal and by (|27|l obey 
[eo,ei] = -6*1161, 
[eo,e2'] = -6*2262', [eo,e3'] = -^2263', (81) 

[ei,e2'] =71.262', [61,63/] ==-71263/, 

[e2',e3'] = 0. (82) 

Following §6 of [43 and noting that n2 7^ because of 
{75, we may choose coordinates i', x', y', z' such that 

eo-|7> ei=n2(0^, 



e'2=We^'|„ 



e'^^B{t')e- 



dz' 



Relations (I81|l and H82|) are satisfied if and only if 



(ilnri2 



d\nB 



— —Oil, 



(83) 



(84) 



dt' ' dt 

By (|77|l . H78|l and (|56|l . the autonomous dynamical sys- 
tem for the fundamental variables On and ^22 becomes 

dOii 



dt' 

de22 
dt' 



^22 — 01l)(6'll + 2^22), 



722 — I'll IC'22- 



(85) 
(86) 



As Oil = 022 and 6*22 = are not allowed by Lemma 1 
(b) and e.g. ^+^, one derives from (jHtj and IJHEI) 
that 



9ii = (21n022- Co + 1)022 or 
u = ^-1^21n022-Co, 

[722 



(87) 



with Co an integration constant. Using u as time vari- 
able, this yields 

022 = ^exp(^), 0ii = ^exp(^)(^ + l), (88) 

with 

C^2exp(^). 

From (|SS|l and l|S7|) one obtains du ~ —2022udt' = 
— Cuexp {^)dt' , such that H84|l may be integrated to give 

n2(u) =Cie^u5, B{u)^C2U^, (89) 

with Ci and C2 non-zero constants. From (|76|l one finds 
8C1 — C^, and after a coordinate change u — 2e~*, x' — 
x/2,y' = ^/2C2y,z' — \/2C2Z the metric reads 

C'^ds^ = exp(-2e"*)(-dt2-|-e*dx2) 

^'e\e-'■'dy^ + e'dz^). (90) 

Thus we find a unique spacetime up to a constant rescal- 
ing. The only non-zero components of the projected Weyl 
and shear tensors are 



i?2'3'(t) = i?22(i) 

= -^exp(--t + 2e-*) 

(T2'2'(i) = 0-3/3/ (i) = --(Til (t) 

= --Cexp(-i-f e"*). 



(91) 



(92) 



while the scalar expansion rate, energy density and pres- 
sure are given by 

0(i) - aexp(e-*)(e-*-)-^), (93) 

Mi) = ^exp(2e-*)(e-* + l), (94) 

Ki) = ^exp(2e-*)(e-*-l). (95) 

Note from (|94|l and (|95(l that constant p is not allowed, 
which could be deduced more directly from H78|l and 
(J55|l - (J5S|) : together with theorem 4.2 this provides an al- 
ternative proof for the non-existence of irrotational PM 
dust [431 . One may eliminate t from (|94|l and (|95|1 via 
[I — p = ^^exp(2e~*) > 0, which yields the following 
equation of state 



M + P=2(A^-P)ln^^^ 



(96) 



Cexp( 



-t\a_ 



which is future directed 



We have u" 

for C > 0. We conclude that (|^ is the metric of a 
perfect fluid model, which starts off with a stiff matter- 
like big-bang singularity at a finite proper time in the 
past (corresponding to t = — oo, with lini(_,_oo ^J'/p =1) 
and which expands indefinitely towards an Einstein space 
with ^i{oo) = ^,p{oo) = -^,9{oo) = ^. Note that 
the dominant energy condition ^ > 0, — /x < p < /j, is sat- 
isfied throughout spacetime, whereas at t — (i.e. after 
a proper time /j°° exp(— u)^) p becomes negative. 
We conclude with: 

Theorem 4.3. Up to a constant rescaling of the met- 
ric, there exists a unique purely magnetic perfect fluid 
which satisfies any two of the three properties of The- 
orem 4.2. This fluid is non-rotating and has vanishing 
spatial gradients of energy density and pressure, with 
the line element given by (|^ and the equation of state 
by (|TO . It is orthogonally spatially homogeneous of 
Bianchi type VIq. The Petrov type is I{M°°) in the ex- 
tended Arianrhod-Mclntosh classification, the shear ten- 
sor being degenerate in the plane perpendicular to the 
0-eigendirection of the projected Weyl tensor. 



4.4. Relaxing the purely magnetic condition. 

We want to indicate here how the algebraically general 
purely magnetic spacetime of theorem 4.2 (and at the 
same Lozanovski's type D class |33) naturally fits into 
a wider class of perfect fluid models. More precisely, we 
drop the purely magnetic condition and look at the class 
A of non- vacuum, non-conformally flat, non-rotating per- 
fect fluids {M,gab,u°') which have a degenerate shear 
tensor aab 7^ and vanishing spatial gradients of energy 
density and pressure. As uja — Ua — DaP = DaH ~ 0, 
one derives DaO = from (|22|l and H12|l with / = /x, 
after which DaicJbcf^^") = follows from ^^ and ipi|) 
with f = 9. The shear being degenerate, abcO'^'^ is the 
only independent scalar which may be built from aab- As 
iia ^ LOa = 0, we may choose a time coordinate t such 
that u = -gr: HtPtO and aabO'"''' depend then on t only and 
it follows that A is part of the class of so called kinemat- 
ically homogeneous perfect fluids, defined and studied in 
[43 . It follows from the analysis there that for any mem- 
ber of A, an eigenframe B = {u,eQ} of aab exists, for 
which H69|l - 171|l and H73|l hold, the shear being degener- 
ate in the (e2,e3) plane {hi — h2 + hs = 0). Herewith, 
it follows from ((TH|l and ((T7|) that aab, Eab and Hab diag- 
onalize simultaneously in B. Thus the fluid congruence 
■u°- is Weyl principal and, as follows from the introduc- 
tion, purely magnetic (purely electric) spacetimes of A 
automatically have a PM (PE) Weyl tensor w.r.t. m". 
From the (22) and (33) components of (|18|l and (|17() one 
obtains 



Also, the difference of the (22) and (33) components of 
the field equations reduces to ni(n2 — n^) = 0. 

For 71,2 = na the Petrov type is D, and it was further 
deduced in 43' that the corresponding models are spa- 
tially homogeneous and LRS III. Note from (|97|l that any 
PE model in A is of Petrov type D (with ^2 ~ n^ ~ 0) 
and was proved ^01 to belong to the Szekeres-Szafron 
family 51, 52, 53, 54]. This family is here characterized 
as the PE subclass of ^, whereas Lozanovski's family |33 
forms precisely the PM type D subclass of A. 

For ni = on the other hand, one reobtains H72(l from 



3 

1 2 



2 
3 1 



2 3 



while Hii = 0, i/33 = -H22 ^ then follows from (j^ZJ. 
For n2 = na = one obtains a subclass of the PE type 
D models of A, while the Petrov type is I if and only 
if n2 7^ 0. Moreover, the essentially unique spacetime 
of theorem 4.3 can now be characterized alternatively as 
the unique type I(M°°) member of A; in this respect, 
also note that the invariants |(SJ) and © become 



J = -2E22{El2+Hl2), M = - 



2H22{9E22 



Hh? 



^-^22(^22 



-^22) 



Further, one derives that qa — r^ — don2 = in 
this case, just as in section ^21 (following 43J). Thus 
the corresponding spacetimes are spatially homogeneous 
(da = 0) of Bianchi type VIq. The surviving equa- 
tions describe the so-called "degenerate shear" subclass 
St{VIo) of the Bianchi Vh family^. 



£^22 = 


2 1 

-nl + -h2{h2 + e), 


(98) 


-^22 — 


^2^2, 


(99) 


M = 


2 S^-hl 

-2+ 3 , 


(100) 


de 

'di ^ 


-li0'' + hl-nl + 3p), 


(101) 


dh2 
dt 


-eh2 - 2nl, 


(102) 


dn2 
dt 


-^{e-2h2)n2. 


(103) 



E2 



E: 



33, 



H2 



-/l2"-3, ^33 = -h2n2. (97) 



The algebraic relations (|98|l - l|100l) are consistent under 
evolution by (|101|l - l|103l) and the Bianchi propagation 
equations. Hence for every choice of the free function 
p{t), there is a family of solutions to the Einstein equa- 
tions corresponding with the dynamical system H101|l - 
(|103|l . When a specific barotropic equation of state 
p = p{fJ-) is assumed, this system becomes autonomous 
by 1)100(1 : the case p = —A is of particular interest since 
a family of irrotational dust spacetimes is then obtained, 
but an explicit integration is not known. Another possi- 
bility to extract a subfamily is to impose conditions on 
the Weyl tensor, e.g. E22 = 6i?22 for constant b. By 
(IHEl and ((^ this equation turns (|102|I - H103|) into an au- 
tonomous dynamical system in /12 and n2, parametrized 
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by b, the exact solution of which can be given in terms of 
elhptic integrals. Alternatively, £'22 = bH22 and its time 
evolution ensure that n2 and p become algebraically de- 
pendent of /i2 and 9, which turns I|1U1|) - (|102|I into an 
autonomous system. The case 6 = eventually yields 
the purely magnetic spacetimes (I90|l . Finally note that 
imposing both a specific equation of state and a rela- 
tion of the form £^22 = 6-^22 , would generically force the 
mentioned dynamical systems to be inconsistent, which 
clarifies again the result for PM irrotational dust jl^ . 



5. CONCLUSION 

Perfect fluid solutions of the Einstein field equations 
were considered, for which the Weyl tensor is purely mag- 
netic with respect to the fiuid velocity u°. We first gen- 
eralized the results |4^ on the non-existence of the so- 
called anti-Newtonian non-rotating dust universes, to the 
case of rotating dust which is either of Petrov type D, or 
which has a vanishing spatial gradient of the matter den- 
sity. These results, in combination with some ongoing 
work on rotating dust (for which the further subcase of 
degenerate shear has by now been dealt with) lead us to 
conjecture that purely magnetic dust spacetimes do not 
exist. 

Motivated by the existence of non-rotating purely mag- 
netic perfect fluids of Petrov type D, for which both 
the spatial gradients of energy density and pressure van- 
ish [33, we studied the algebraically general case and 
demonstrated that any two of the following conditions 
implies the third: (i) the fluid congruence is non-rotating, 
(ii) the fluid is non-accelerating, (iii) the fluid has a van- 
ishing spatial gradient of the matter density. For purely 
magnetic perfect fluids satisfying these conditions it turns 
out that the magnetic part of the Weyl tensor has one 
eigenvalue and that the shear tensor is degenerate in the 
plane orthogonal to the corresponding eigenvector. The 
unique (up to constant rescaling) solution satisfying these 
conditions is an orthogonally spatially homogeneous per- 
fect fluid of Bianchi type VIq which has a big bang sin- 
gularity, starts of as stiff fluid and asymptotically evolves 
towards an Einstein space. 

Finally notice that purely magnetic spacetimes of 
Petrov type D or I{M°°), as considered in the present 
work, cannot be conformally mapped to vacuum space- 
times, by the results of Q and [ll| and the fact that 
the Weyl tensor is preserved under such mappings. The 
question remains open whether type I{M+) purely 
magnetic spacetimes exist at all. 

An overview of the results obtained in this paper for 
Petrov type I is given in table I. 

Note added in proof: it has been brought to our 
attention that the same metric of section 4.3 has been 
independently found by C. Lozanovski. 
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APPENDIX 

We provide here the quite long and technical proof 
for the part (i), (ii) => (Hi) of proposition 4.2, stat- 
ing that for non-rotating and non-accelerating type I 
PMpf 's the spatial 3-gradient of the energy density van- 
ishes. Projection w.r.t. an orthonormal eigenframe B of 
aab turns out to be favorable, and we actually prove 
the result regardless of the Petrov type. We will split 
the proof in two cases: degenerate and non-degenerate 
shear. For conciseness we will write p = fi + p and we 
add TTia = dafJ. = daP, Zq = daO as variables. Note that 
we may exclude p = by the result of [Tjl or [13 . 

With CTq+iq-i — 0, the (a-|-la— 1) component of l(TH|) 
reads Waha — 0, such that Wq = when the shear tensor 
is non-degenerate; for degenerate shear, say hi = 0, only 
W2 = W3 = follows, but we may then partially fix 
the frame by taking Wi = (hereby leaving the freedom 
of performing rotations about an angle x in the (62, 63)- 
plane satisfying Oqx — 0), such that also here the triad 
{e„} is Fermi-propagated along the fluid flow |43- Hence 
the equations H53l) - (|61|l are valid. Further, the Bianchi 
constraint divE^ translates to 



'ShaHa^ 



-1=0. 



(A.l) 



On applying do to Qa and eliminating Ha+ia-i by 
means of Qa , one gets 



donia 



5 1 

Om„ + -{ha+i - ha-i)mo 



(A.2) 



3 "6 
Herewith, [do,da]p hecovaes 

Ra^ pza~ -{26 + ha + i- ha-l)ma=0 (A. 3) 

6 

while [do, da]0 results in an expression for doZa'. 

1 



doz, 



'QZa 



-9Za + {ha + 1 — ha-l)Za + 
-ZQah^^l — 2ro,ha^i 



6 



(A.4) 

(I) Suppose the shear tensor is degenerate, say hi — 
0. Then /ig = -ft.2 7^ 0, mi = by i?i, while 
{Q2,R2,d2hi = 0} (resp. {Q^^Rs^dsihi) = 0}) can 
be solved for Hi^, 22, 92 (resp. H12, Z3, r^) to give: 



^13 — 



92 



1712 

~3^' 
m2 (2/126* 



Z2 



m2{2e-h2) 



6p 



hi 



2p) 



6phl 



H12 



rs 



ma 
3/12' 



Z3 



m3{29~h2) 
Qp 



(A.5) 



(A.6) 



7713(2/12^ -hl + 2p) 
Qphl 
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p = const p 7^ const Dfi 7^ Dfi — 
A $ conjecture: $ 

iP|l $ 

TABLE I: Results obtained for PMpf's of Petrov type I: $ and A indicate respectively that no solutions exist or that all 
solutions are known. 



Computing the time evolution of R2, resp. R3 and in- 
serting HA.5|I . resp. HA.6|) . one gets 



1712 Oi — ni3 Oi = 



with 



Oi = -2p/i2 + 4^6* + 12p2 + 4ph20 



phi 



(A.7) 



4pp. (A.8) 



Suppose now 7712 =/= 0, then Oi = 0. Applying twice 
82 to this equation, in each step making use of l|55(l 
and d2P = d2P — (the latter following from [9o, da\p, 
substituting HA.6(I and dividing by TO2, one gets a linear 
system Oi = O2 = O3 = in the two variables p,p, 
parametrized by p, 0, /12; the vanishing of the extended 
system determinant yields 



Pi = -9/i^ + b4ehl + 274hlp - 108r/i2 - 648/i^p6' 



-726'^/ii 



168/i^6i2p- 



-5526ip2/i2 ~ 360p3 



h 660p2/i2 + Mpe^h2 
192p^9^ = 0. 



Repeating the above process for Pi instead of Oi , one 
obtains two new polynomial equations P2,-p3- Elim- 
inating 9 between Pi and P2, resp. Pi and P3 by 
calculating the resultant w.r.t. 9 yields two equations 
hlp^Ki{p,9) = and hlp^K2{p,9) = 0, where the re- 
sultant of Ki and K2 w.r.t. p is proportional to ftl*^. 
Hence one obtains the contradiction p = or /12 = 
(lemma 4.1 (b)). Thus m2 = 0. The same reasoning, 
involving the same polynomials, holds when applying 
^3 instead of 82 in the above, such that also 7713 = 0. 
Hence for degenerate shear, the fluid has a vanishing 
spatial gradient of the energy density. 

(H) Suppose the shear is non-degenerate, i.e. /11/12/13 7^ 
0. On making use of (|54|l and H55|l one 

can solve {[92, 93]/9, [^2, 93]6', 92Q3, 93(92} for 
92TO3, 93W2, i922;3, 93^2, yielding 



927713 = (737712 - -r-niimi 
hi 

3 1 

2p 3/l2^3 



7ri2TO3. 



(A.9) 



Then [9o,92]7n3 implies 

Xi = 9hih2h3m2m3p -\ 
+3h2h^niimi — ft.i(/i2 



[3/l2^3/il('"2"73 + 937772) 
- /73)77l27773]p^ = (A.IO) 



where in (jA.Qp and ljA.10|) we have eliminated the 
Ha+i,a-i by Qa tcrms and the Za by Ra- By repeat- 
ing the above for the index couples 31 and 12 instead 
of 23, one gets equations ^2,^3, which can formally 
be obtained by cyclic permutation of the indices of Xi 
(twice). Taking the combination 7772 A'2 — 7773X3 one 
gets (with p ^ 0): 



Aq = 3hih^m2n22 — 3ft.ift.2'™3'^33 

+ 3/ll/l2/l3('7727773(/l2gi " hsVl) 
+777l7772/l2?'3 - "'-l?n-3/''392) 
= 777l7772n73(/li -f 2/12/13). 



(A.ll) 



On examination of (|55|l . l(5^ . l(HT|) and IIA.211 one sees 
that calculation of A^+i = doAi + - — ^-^9Ai,i = 
0,...,5, yields a system {Ai,i — 0..6} of 7 hnear 
equations in 7722, 7733, qi,ri,q2,r^, parametrized by the 
ma and /ic, which can only be consistent if the de- 
terminant of the extended system matrix vanishes. 
This yields mlmlmlh\^hfhfCi{h2,h-i) = 0, where 
Ci(ft.2,/73) is a homogeneous polynomial of degree 8; 
thus if Ci{h2,h^) vanished the ratio ^12/^3 would 
be constant, in contradiction with lemma 1 (c) and 
the assumed non-degeneracy of the shear. There- 
fore 771177727773 = 0. As a vanishing spatial gradient 
of the energy density (all 777i zero) implies degener- 
ate shear (cf. section 4), we are left with two qual- 
itatively different subcases: (a) 7771 — 0,77727773 7^ 
and (b) 7772 = 7773 = 0, 7771 7^ 0. In general, if 
nia = for fixed a, it follows that Ha+ia-i = Za ^ 
from (|A.1|I and HA.3I) . From (|A.4|) . (|5S1) and (jHU one 
then deduces La = /i^_|_i9a — h^_ira — 0, doLa = 
{4:ha+i + 5ha-i)hl_^_iqa + {4:ha-i + 5ha+i)hl_ira = 0, 
which for non-degenerate shear can only be consistent 
if qa = ra= 0. 
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(a) mi — 0, 17121123 ^ 0: 

As also Qi = ri = 0, one deduces 7122 = ^33 = 
from X2 = X3 = 0, i.e., 112 — n^ — —iii = -^''^n. 
Then e.g. the (31) component of the field equa- 
tions together with (-190) ^^^ the (31) com- 
ponent of (|18|) together with the (2) component 
of (|24(l respectively give 

92^1 == -2g2?^l, 
52^^33 = iH33~H22)r2-SHuni. (A.12) 



while (I53|) further reduces to 

H22 =nih2, H33 = -nihj 



(A.13) 



Propagation of the second equation of 
along the 82 integral curves, using if^ . 
and HA.13I) . yields ?ii(4i/i3 + Z2 + 4/13^2) = 0. 
As 7ii = is not allowed by ljA.13|l . it follows 
that Bq = 4i/i3 + Z2 + 4/13^2 — 0. Calculation of 
Bi+i = doBi + ^^OBi, i = 0, 1, 2 gives a homo- 
geneous Hnear system {Q2, -Bq, 51,-62,-83} in the 
variables 92, ''27 '^^^2, -^2, -^'^13 and parametrized by 
/i2,/i3, which can only be consistent with TO2 7^ 
if the determinant, namely |ft,^/i2^§(21/i2 + 
64/13/12 + 16/13), vanishes, which is contradictory 
to lemma 4.1 (c). 

(b) 7712 — 1713 = 0, mi 7^ 0: 

Apart from q2 = ^2 = Qs = r^ ^ 0, one has now 
Till — from (jA.lOp . One can solve the (23) com- 
ponent of the field equations together with its do 
derivative and [do, i9i 1(^22 - "33) for i9in22, 9in33 
to give 

c'i"-22 = -777— [(4i?23 - zi+ Ah2ri)n22 
Ztli 

+ {4H23 + zi+ 4h2qi)n33] , (A.14) 

9in33 = -T-T— [(4if23 - ^^1 -4/l3''l)"-22 
Ztli 

+ {AH23 + zi - 4/1391)7133] . (A.15) 



Then [5o,9i](n22 + "33), eliminating i?23 by i?i, 
gives 



Co = 12 [/l2n22 - (2/l3 -f- 3/l2)"-33] ^291 

-1-12 [/l3"-33 - (2/l2 + 3/13)7122] /l3''l 
-1-13(7122 - "33)™l - 6 [(/l2 + 3/13)^22 

+ (/i3 + 3/72)^33] ^1 = 0. (A.16) 



Calculation of Ci- 



doa 



6+2i 



oa 



0, ... ,3, eliminating in each step iJ23 by means 
of Ri, gives two homogeneous linear systems 
{Co,Ci,C2,Cj}, j — 3,4, in the variables 
gi, r I, Till, zi and parametrized by /12, /73, 7122, '^•33- 
Again this can only be consistent with 7771 7^ 
if the determinants of the coefficient matri- 
ces vanish. However, on computing the re- 
sultant of these determinants w.r.t. 7722, resp. 
7133, one gets equations hf'h\h^P{h2,h3)n^ — 
/if /i|/i|P(/i2, h3)nll = 0, with P(/72, /13) a homo- 
geneous polynomial of its arguments, such that 
7i22 = "-33 = by lemma 4.1 (c). Hence 771 — 
71-2 = 713 — 0, in contradiction with (|53|) . 



Hence case (II) is not allowed, and this concludes the 
proof. 



[1] A. Matte Canadian J. Math. 5, 1 (1953) 

[2] L. Bel Cah. Phys. 16, 59 (1962) (English translation Gen. 

Rel Grav. 32, 2047 (2000)) 
[3] J.J. Ferrando and J. A. Saez Class. Quantum Grav. 19, 

2437 (2002) 
[4] R. Maartens and B.A. Bassett Class. Quantum Grav. 15, 

705 (1998) 
[5] A. Barnes Proceedings 27th Spanish Relativity Meeting, 

Alicante, Spain. Sept. 2003, eds. J. A. Miralles, J. A. Font 

and J. A. Pons (Univ. Alicante Press, 2004), see also gr- 

qc/0401068 
[6] C. Lozanovski and J. Carminati Gen. Rel. Grav. 34, 853 

(2002) 
[7] C.B.G. Mcintosh and R. Arianrhod Class. Quantum 



[9] 
[10] 

[11] 
[12] 

[13] 
[14] 
[15] 
[16] 



Grav. 7, L213 (1990) 

R. Penrose and W. Rindler Spinors and Spacetime vol 2 

(Cambridge, Cambridge University Press, 1986) 

C.B.G. Mcintosh, R. Arianrhod, S.T. Wade and C. 

Hoenselaers Class. Quantum Grav. 11, 1555 (1994) 

R. Arianrhod and C.B.G. Mcintosh Class. Quantum 

Grav. 9, 1969 (1992) 

C.H. Brans J. Math. Phys. 16, 1008 (1975) 

B.M. Haddow J. Math. Phys. 36, 5848 (1995) 

M. Triimper J. Math. Phys. 6, 584 (1965) 

N. Van den Bergh Class. Quantum Grav. 20, LI (2003) 

N. Van den Bergh Class. Quantum Grav. 20, L165 (2003) 

J.J. Ferrando and J. A. Saez Gen. Rel. Grav. 36, 2497 

(2004) 



13 



[it; 
[is; 

[19 

[20 

[21 

[22; 

[23; 

[24' 
[25^ 
[26; 

[27; 

[28; 

[29 

[30 
[31 

[32; 

[33 

[34; 

[35; 

[36 

[37' 
[38; 

[39 

[40 

[41 
[42; 



Carminati Gen. Rel. Grav. 37, 


605 


[43] 
[44] 


.A. Saez Class. Quantum Grav. 


20, 


[45] 


J. A. Saez J. Math. Phys. 45, 


652 


[46] 



E. Zakhary and J. 
(2005) 

J.J. Ferrando and J. 
2835 (2003) 
J.J. Ferrando and 
(2004) 

R. Arianrhod, A.W-C. Lun, C.B.G. Mcintosh and Z. [47] 
Perjes Class. Quantum Grav. 11, 2331 (1994) 
R. Maartens, W.M. Lesame and G.F.R. Ellis Class. 
Quantum Grav. 15, 1005 (1998) [48] 

C. Lozanovski and J. Carminati Class. Quantum Grav. 
20, 215 (2003) 

J.M. Stewart and G.F.R. Ellis, J. Math. Phys. 9, 1072 [49] 
(1968) 

S.W. Hawking Astrophys. J. 145, 544 (1966) [50] 

G.F.R. Ellis and M. Bruni Phys Rev D 40, 1804 (1989) 
P.K.S. Dunsby, B.A. Bassett and G.F.R. Ellis Class. [51] 

Quantum Grav. 14, 1215 (1997) [52] 

R. Maartens, G.F.R. Ellis and S. Siklos Class. Quantum [53] 
Grav. 14, 1927 (1997) 

G.F.R. Ellis and P.K.S. Dunsby Astrophys. J. 479, 97 [54] 
(1997) 

C. Lozanovski and M. Aarons Class. Quantum Grav. 16, [55] 
4075 (1999) 

W.B. Bonnor Class. Quantum Grav. 12, 1483 (1995) 
C. Hillman Electrogravitism versus Magnetogravitism [56] 
|http://math . ucr.edu/home/baez/PUB/electromagneto 
C. Lozanovski and C.B.G. Mcintosh Gen. Rel. Grav. 31, 
1355 (1999) 

C.B. Collins and J.M. Stewart Mon. Not. R. Astron. Soc. 
153, 419 (1971) 

G. Fodor, M. Marklund and Z. Perjes Class. Quantum 
Grav. 16, 453 (1999) 

H. van Elst and G.F.R. Ellis Class. Quantum Grav. 13, 
1099 (1996) [57] 

M. Marklund Class. Quantum Grav. 14, 1267 (1997) 
C. Lozanovski Class. Quantum Grav. 19, 6377 (2002) [53] 

N. Van den Bergh and L. Wylleman Class. Quantum 
Grav. 23, 3353 (2006) 

S. Matarrese, O. Pantano and D. Saez Phys. Rev. Lett. 
72, 320 (1994) 

H. van Elst, C. Uggla, W.M. Lesame, G.F.R. Ellis and 
R. Maartens Class. Quantum Grav. 14, 1151 (1997) 
C.F. Sopuerta Phys. Rev. D 55, 5936 (1997) 
L. Wylleman Class. Quantum Grav. 23, 2727 (2006) 



N. Van den Bergh Class. Quantum Grav. 6, 797 (1989) 
R.M. Wald General Relativity{The University of Chicago 
Press, 1984) 

R. Maartens Phys. Rev. D 55, 463 (1997) 
G.F.R. Ellis General Relativity and Gosmology, edited by 
R. K. Sachs (Academic, New York, 1971) 
D. Kramer, H. Stephani, M.A.H. MacCallum, C. Hoense- 
laers, E. Herlt, Exact solutions of Einstein's field equa- 
tions (Cambridge, Cambridge University Press, 2003) 
M.A.H. MacCalluraCosmological Models from a Geomet- 
ric Point of View (Gargese) Vol 6 p 61 (New York, Gor- 
don and Breach, 1971) 

G.F.R. Ellis and M.A.H. MacCallum Commun. Math. 
Phys. 12, 108 (1969) 

A. Barnes and R. Rowlingson Class. Quantum Grav. 6, 
949 (1989) 

P. Szekeres Comm. Math. Phys. 41, 56 (1975) 
D.A. Szafron J. Math. Phys. 18, 1673 (1977) 
D.A. Szafron and J. Wainwright J. Math. Phys. 18, 1688 
(1977) 

J.D. Barrow and J. Stein-Schabes Phys. Lett. 103A, 315 
(1984) 

J. Wainwright and G.F.R. Ellis Dynamical Systems in 
Cosmology chapter 6 (Cambridge, Cambridge University 
Press, 1997) 

However, for LRS class I a complicated third order ODE 
needs to be solved. For LRS class HI the metric was 
shown to be explicitly constructable, but a simpler form 
than the one suggested in [23I can be obtained by di- 
rectly taking the function s(t) as a new time-variable 



u instead of t\ the result is C ds 



i f-idu^ +u^x 



[dx 4- qe^{zdy - 



■ydz)]^ + e^^{dy^-\-dz^)), where C,c 



and q are constants and k equals -1, or 1, i? = 
- log (1 -\- 1^2 _^ ^2^^ ^^^ p ^ 2q'^^4 _ ^^2 ^ ^ 

Following [i^. different from the convention in J47| . 
where J71230 = — 1 is taken. 

The ila are the non-zero components of the angular ve- 
locity vector of the triad {bq} w.r.t. the 'inertial com- 
pass', see |4q| and references therein. The notation Uaa 
is also in agreement with [i^. where the further decom- 
position 7"^^ = e/5c«n«" -H S^an - 5'jla^, n°''^ = n'"''' of 
the purely spatial coefBcients is exploited. 



